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While the stoichiometric intermetallic compound FeAl is found to be paramagnetic in experiment, 
standard band-theory approaches predict the material to be ferromagnetic. We show that this 
discrepancy can be overcome by a better treatment of electronic correlations with density functional 
plus dynamical mean field theory. Our results show no ferromagnetism down to 100 K and since the 
susceptibility is decreasing at the lowest temperatures studied we also do not expect ferromagnetism 
at even lower temperatures. This behavior is found to originate from temporal quantum fluctuations 
that screen short-lived local magnetic moments of 1.6 /tb on Fe. 

PACS numbers: 71.10.-w, 71.27.-|-a 


I. INTRODUCTION 

Intermetallic alloys of iron and aluminum have a high 
hardness with a much lower specific weight than steel. 
Because of this, their low costs, and resistance against 
corrosion and oxidation, FeAl alloys are often used as 
lightweight structural materials. Most puzzling are the 
magnetic properties. Here, experiments such as high- 
held Mossbauer investigations^ indicate no magnetism 
for stoichiometric FeAl which forms a B2 GsCl-type of 
lattice (two interpenetrating Fe and A1 simple cubic lat¬ 
tices) . Especially the fact that FeAl does not show ferro¬ 
magnetism in experiment, while electronic structure cal¬ 
culations within spin-polarized density functional theory 
(DFT) predict a ferromagnetic ground state has drawn 
attention to the material: independently of the band- 
structure code, DFT orbital basis set and exchange corre¬ 
lation potential a ferromagnetic ground state with a mag¬ 
netic moment at the Fe site of about 0.7 is found. 
Even though the energy difference between the ferromag¬ 
netic and the nonmagnetic state is rather small, the fer¬ 
romagnetic state is stable over a wide volume range. In 
fact, only a reduction of the lattice constant by more than 
10% would suppress ferromagnetism.® This high stability 
of the ferromagnetic phase in FeAl suggests that the devi¬ 
ation from experiment is not just a numerical inaccuracy, 
but requires a deeper understanding. 

Different approaches have been used hitherto to ex¬ 
plain the deviation between spin-polarized DFT and ex¬ 
periment. One explanation is based on the fact that the 
processes used to prepare FeAl often ’’freeze in” chem¬ 
ical disorder. That is, ’’real” FeAl is usually not fully 
ordered due to various lattice defects, such as vacancies 
and antisites, which in turn could have a significant effect 
on the magnetic properties of the material. Against this 
background, there exist several studies concerning the ef¬ 
fects of disorder on the magnetic properties of FeAl.®’®’^’® 
For example, in Ref. 7 the disorder is included via the 
coherent potential approximation (CPA)® in the Kor- 
ringa, Kohn and Rostoker (KKR) framework,with 
the paramagnetic phase described by the disordered lo¬ 


cal moment approximation (DLM).^® In agreement with 
previous DFT calculations, it has been found that ideal 
FeAl is ferromagnetic. However, even with a small degree 
of disorder the paramagnetic state, without net magne¬ 
tization but nonzero local moments, becomes the stable 
configuration. Thus, disorder destroys the long-range fer¬ 
romagnetic order in DFT. 

However, no ferromagnetism has ever been observed 
for stoichiometric “real” FeAl, even for samples with 
very low defect concentration. Therefore, it still remains 
the question if perfectly ordered FeAl would really be 
ferromagnetic as predicted by DFT. Indeed, Mossbauer 
experiments^ find magnetic moments only for Fe anti¬ 
structure atoms (which means Fe atoms sitting on an A1 
lattice site) and their eight Fe neighbors. 

Another possible explanation has been given in Ref. 6 
using the DFT-I-C/ approach.^® Usually, one would expect 
DFT-I-I/ to yield larger magnetic moments and a stronger 
tendency towards ferromagnetism than DFT. For U val¬ 
ues ranging from 4 to 5 eV a nonmagnetic state however 
coexists with the ferromagnetic one in DFT-I-I/. The 
ferromagnetic state even disappears for a rather large 
U = 5eV, which offers another explanation of the non¬ 
magnetic nature of FeAl. This rather unusual DFT-I-I/ 
result can be explained by the changes in the density of 
states (DOS): increasing U reduces the DOS at the Fermi 
level so that according to the Stoner criterion there is no 
ferromagnetism® even though the effective exchange is 
increased by U Hence, in a narrow range of U, there 
is no ferromagnetism in DFT-l-I/.® 

In Ref. 14 it has been argued that this DFT-I-I/ result 
has to be taken with a grain of salt and it has been pro¬ 
posed for the first time that dynamical spin fluctuations 
suppress ferromagnetism in FeAl. This has been sup¬ 
ported by a dynamical mean field theory (DMFT^®’^®) 
calculation.^^ For U = 2eV, FeAl is found^^ to be para¬ 
magnetic in DFT-|-DMFT.^^“®® However, Ref. 14 only 
shows a single DFT-I-DMFT result, the spectral function. 
The proposed spin fluctuations, the magnetic properties 
and susceptibility have not been calculated. 

Considering these limited results as well as the im- 
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provements of DFT+DMFT in recent years, a more thor¬ 
ough analysis is in order. Beyond the first DFT-I-DMFT 
spectrum of Ref. 14, we study the local and bulk magnetic 
susceptibility, the magnetic moment and the fc-resolved 
spectrum. We also explicitly calculate the local inter¬ 
actions ab-initio by constrained random phase approxi¬ 
mation (cRPA) and beyond Ref. 14 we include the cal¬ 
culated Hund’s exchange in DMFT with its full SU(2) 
symmetry, since it plays a pivotal role for the magnetic 
properties. Our results show that while there is a local 
moment of even 1.6 /tb on short time scales, it is screened 
(suppressed) on longer time scales. This suppression of 
the local moment occurs on the fs time scale (eV“^) and 
explains why there is eventually no long-range ferromag¬ 
netic order. 


In Section II we present the DFT bandstructure and 
DOS as well as the Wannier function projection. Section 
III is devoted to the one-particle properties as calculated 
in DFT-I-DMFT, i.e., the self-energy as well as the lo¬ 
cal and fc-resolved spectral function. The DFT-f-DMFT 
magnetic properties are discussed in Section IV, i.e., the 
local and (zero) ferromagnetic moment as well as the 
time-dependent local susceptibility and bulk susceptibil¬ 
ity. Finally, Section V summarizes the results and puts 
them into context with experiment. 


II. ELECTRONIC STRUCTURE WITHIN DFT 


As a first step, we employ the Vienna ab initio Sim¬ 
ulation Package (VASP)^^ with GGA-PBE functional^^ 
for calculating the bandstructure and density of states of 
FeAl. Fig. 1 shows the bandstructure of FeAl around 
the Fermi level. The bands closest to the Fermi level 
have mainly Fe 3d character and are split into t 2 g and 
Cg due to the cubic crystal field. For these bands we 
will later include electronic correlations by DMFT. How¬ 
ever, since the Fe 3d bands strongly hybridize with the 
A1 3s and 3p states, we also include these A1 bands (as 
non-interacting) in our low energy Hamiltonian. The cor¬ 
responding Hamiltonian is obtained by a projection onto 
nine maximally localized Wannier orbitals^^, which re¬ 
produce the DFT bandstructure well, see Fig. 1. 



FIG. 1: Electronic bandstructure of FeAl: Gray points denote 
the DFT bandstructure while colored lines show that of the 
Wannier projection. The color of the bands indicates the 
amount of Fe t 2 g (red), Fe Cg (green) and Al-sp^ (blue) orbital 
contribution. 
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FIG. 2: FeAl density of states obtained by DFT 

Fig. 2 shows the orbital resolved density of states. 
It can be seen that the central peak is mainly of t 2 g 
character. As the Fe t 2 g orbitals form only weak bonds 
with the nearest-neighbor A1 atoms, they have a rather 
small energy dispersion. The Fe eg states instead point 
towards the neighboring Fe atoms and hybridize more 
strongly. Hence they have a larger bandwidth and split 
into a bonding- and an antibonding-like part. 


III. DMFT SELF ENERGY AND SPECTRAL 
FUNCTION 

After obtaining the low energy Hamiltonian in the ba¬ 
sis of Wannier functions, we perform DMFT calculations 
including the five Fe d orbitals and the 4 A1 sp^ orbitals 
within a so-called dp model.We supplement the DFT- 
based Wannier Hamiltonian in DMFT by a local d-d 
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Kanamori interaction, but disregard d-p and p-p inter¬ 
actions beyond what is already contained in DFT. Note 
that the hopping terms of the Hamiltonian still contain 
the full information about the hybridization with the A1 
sp^ states and charge transfer between d and sp^ orbitals 
is allowed. 

We calculate the screened many-body Coulomb inter¬ 
actions [/, U' and J by the constrained random phase 
approximation (cRPA),^®’^® where we exclude only the 
Fe d states from the screening. This is appropriate as 
interactions are also applied only to these d states. 
For our DMFT calculation, we use the average values 
for the intra-orbital Coulomb interaction U = 3.36 eV, 
the inter-orbital Coulomb interaction U' = 2.36 eV and 
the Hund’s coupling J = 0.71 eV. This yields a local, 
5'?7(2)-symmetric Kanamori interaction;^®’^® 


H\qc — ^ ^ U 

m 

+ E iu' ,—<T (U a] 

m^m' ,(T 

m^m' 

+ E (1) 

m^m' 


Here, (cma-) creates (annihilates) an electron with 
spin fj in the Fe 3d orbital m; rima = clna^ma- We em¬ 
ploy the double counting correction of the fully localized 
limit,®® and validate that a difference of 2.5 eV in the dou¬ 
ble counting does not change our findings (not shown). 

For the solution of the DMFT impurity problem we 
use a continuous-time quantum Monte Carlo (CT-QMC) 
algorithm in its hybridization expansion (CT-HYB) in 
the version of Ref. 31, for a review see Ref. 32. Espe¬ 
cially with regard to the magnetic properties that we will 
compute, it is important to employ the rotationally in¬ 
variant form of the interaction term iJioc above, including 
a pair-hopping and a spin-flip term, and not only density- 
density contributions. As for the CT-HYB, we note that 
it is essential to truncate the outer states for the evalu¬ 
ation of the local fermionic trace only at high energies, 
especially at high temperatures. 

Fig. 3 shows the imaginary part of the DMFT self en¬ 
ergy on the Matsubara axis for all five Fe d orbitals. 
In order to avoid all uncertainties related to an analyt¬ 
ical continuation, we calculate the quasiparticle weight 
Z directly from the self energy on the Matsubara axis 
Z = 1/(1 — 3(9I](*u;)/9(*a;))|itd->-o)- This yields a value 
oi Z = 0.75 , essentially the same for all 3d orbitals. This 
Z value would indicate a rather weakly correlated mate¬ 
rial. 
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FIG. 3: DMFT self energies for the Fe d orbitals (at inverse 
temperature d = 30eV“^ corresponding to 390K). The ex¬ 
tracted quasiparticle weight is Z = 0.75 . 
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FIG. 4: DMFT spectral function for the Fe d orbitals at d = 
30eV“^, compared to DFT (d^; = 0 corresponds to the Fermi 
level). 


The corresponding spectral function A(k, w) = 
— l/7r3(G'(k, w -h fO®')) is shown in Fig.4 on the real fre¬ 
quency axis, for which an analytic continuation using a 
stochastic version of the maximum entropy method has 
been used.®® In comparison to the DFT DOS, both oc¬ 
cupied and empty states are slightly shifted towards the 
Fermi energy due to the Fermi-liquid renormalization. 
There is no evidence of pronounced upper and lower Hub¬ 
bard bands and one can only observe a weak increase of 
the spectral weight at high frequencies. In agreement 
with Ref. 14, we find that the spectral function at the 
Fermi level is essentially the same in DFT-l-DMFT as 
in DFT. We did not perform charge-self-consistent cal¬ 
culations since the difference in the occupation of the 
d-orbitals between the DFT-derived Hamiltonian and 
DMFT is very small. In DFT, we have 4.8 electrons 
in the t 2 g and 2.5 electrons in the Cg states out of 11 
electrons per unit cell, in DMFT the t 2 g orbitals are oc- 



















4 


cupied with 4.8 and the Cg orbitals with 2.6 electrons. 
Also the changes in the one-particle spectrum are rather 
small. Note, only if DMFT alters the spatial charge dis¬ 
tribution p(r) considerably, charge self-consistency would 
have an effect. Thus, we expect changes by charge-self- 
consistency to be small. 

Fig. 5 presents the corresponding A:-resolved spectrum 
which shows that also the DFT-I-DMFT bands essentially 
follow the DFT bandstructure. The most noteworthy 
effects are again a slight shift towards the Fermi level, 
i.e., a quasiparticle renormalization and a broadening of 
the bands, especially of the d bands located around the 
Fermi level. Hence, regarding only single-particle quanti¬ 
ties, FeAl seems to exhibit only weak correlation effects. 
However, this picture changes when considering also two- 
particle quantities, namely the magnetic susceptibility. 



FIG. 5: DMFT fc-resolved spectral function (blue) compared 
to the DFT bandstructure (black). 


IV. DMFT MAGNETIC PROPERTIES 

In order to study the magnetic properties of FeAl 
within DFT-I-DMFT, we compute the local magnetic sus¬ 
ceptibility, represented by the two-particle spin-spin cor¬ 
relation function 

X1oc(t) = xZc^ir) = E iSTir)S^m (2) 

m,n m^n 

with m and n being the orbital indices of the five 
Fe d orbitals, r the imaginary time, and g k, 2 the 
gyromagnetic factor for the electronic spin. = 

l/2{nmtiT) ~ ^m 4 ,('r)) is the z-component of the spin op¬ 
erator of orbital m, expressed in terms of the correspond¬ 
ing density operators rima = cZa-Cma- 

Technically speaking, xioc('r) is obtained by first 
measuring the generalized magnetic susceptibility 
X\oc(.Wtw', iu}) of the converged DMFT impurity model 
by means of CT-HYB quantum Monte Carlo sampling. 
Thereby, y;ioc(*;z, ii/', feu) automatically contains all 


vertex corrections to the bare (DMFT) bubble spin 
susceptibility. The sum over the fermionic Matsub- 
ara frequencies v and v' and a Fourier transform, 
X 1 oc(t) = e-“'"x(*w), finally lead to x 1 oc(t). 

Here, for the large frequency asymptotics, the bare 
bubble contribution, Eq. (3), which is known on a larger 
frequency grid and an additional fitting function of the 
form l/i/^ have been used. The results for xioc(t) are 
shown in Fig. 6 for (3 = 30eV“^ (for lower temperatures 
the numerically feasible frequency box becomes too 
small). 



012345678 

x[eV‘''] 

FIG. 6: Local magnetic susceptibility yioc as a function of 
(imaginary) time r for /3 = 30eV“^. Besides the total sus¬ 
ceptibility, also its diagonal and bare-bubble contribution are 
shown, as well as (in the inset) its Cg and t 2 g contribution.'^® 

The solid, blue curve in Fig. 6 corresponds to the total 
magnetic susceptibility Xioc('r) of Eq. (2). The dashed, 
orange curve instead represents the orbital-diagonal con¬ 
tribution ('’’)■ dotted, purple curve is the 

bare-bubble contribution X?oc(''’)i "'^fiich neglects vertex 
corrections and is obtained by directly convoluting the 
DMFT Green functions Gmiiv)'. 

Xloc(*^) “ 'y ^ 3- iuj) (3) 

iu^ma 

The significant difference between the bare-bubble con¬ 
tribution and the susceptibility including vertex correc¬ 
tions in Fig. 6 reveals that electronic correlations actu¬ 
ally play a major role in FeAl, more than it could be ex¬ 
pected from the weak quasiparticle renormalization. Fig. 
6 also shows that the enhancement of x 1 oc(t) stems ap¬ 
proximately in equal parts from an enhancement of the 
intra-orbital contribution (diagonal part) and additional 
inter-orbital (off-diagonal) contributions, which are not 
present in the bare-bubble susceptibility. 

The local susceptibilities in Fig. 6 show a rather fast 
and strong decay in r. Here, the value of x 1 oc(t) at t = 0 
can be interpreted in terms of the instantaneous, local 
magnetic moment. The observed decay in r reflects a 
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dynamical screening of this local magnetic moment due 
to quantum fluctuations. Thus, we can conclude that 
dynamical quantum fluctuations significantly reduce the 
local magnetic moment in FeAl. 


Fitting xioc('r) to an exponential between r = OeV“^ 
and r = 5eV“^ yields a time scale for the screening of 
Tg = 1.03 eV“^ = 4.02 fs. The inverse of Tg is the energy 
scale associated with the screening which is essentially 
the bare bandwidth if we have a noninteracting system, 
the width of the central peak if we consider the inter¬ 
acting bubble, and the Kondo temperature for the inter¬ 
acting system with vertex corrections. This Kondo tem¬ 
perature is smaller than the width of the central peak."'^ 
Hence the decay with vertex corrections should be slower. 
Indeed, in Fig. 6 the total xioc(t) decays slower than the 
bubble contribution. For a related analysis, how to in¬ 
terpret the susceptibility as a function of imaginary time 
and how the local, fluctuating magnetic moment reflects 
as a pronounced low-energy peak in the local neutron 
spectra, see 34-36. 
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FIG. 7: DFT-(-DMFT magnetization for different tempera¬ 
tures. The ferromagnetic moment is zero within the error 
bars, FeAl is a paramagnet. 


In the inset of Fig. 6, we separate the Cg and t 2 g con¬ 
tributions of the susceptibility. These two contributions 
are rather independent as one clearly sees from the longer 
time scale on which the t 2 g susceptibility decays. This 
different decay rate can be explained by the considerably 
more narrow t 2 g bandwidth and hence stronger corre¬ 
lations of the t 2 g orbitals. If Hund’s exchange was the 
major player, on the other hand, one would expect a 
stronger coupling of Cg and t 2 g susceptibility, and a de¬ 
cay on a similar time scale. 


This all suggests that the Hund’s rule exchange J, 
which mainly drives the inter-orbital contribution, is not 
exceedingly important in FeAl. This is in contrast to 
other Fe-based compounds such as the iron-based super¬ 
conductor LaFeAsO.^^’^^ which have been classified as 
Hund’s metals. 


From the local magnetic properties, we now turn to the 
bulk magnetic susceptibility and the long-range ordered 
ferromagnetic moment. Fig. 7 shows the ordered mag¬ 
netic moment, which has been obtained by breaking the 
spin symmetry in the first DMFT iteration so that the 
system can either stabilize a para- or ferromagnetic solu¬ 
tion. As Fig. 7 clearly shows, the ordered ferromagnetic 
moment is zero down to a temperature of 100 K. Thus, 
in the investigated temperature range, FeAl is paramag¬ 
netic in DFT-I-DMFT, in agreement with experiment but 
in contrast to DFT. 


This result is also supported by the calculation of 
the bulk ferromagnetic susceptibility in DFT-I-DMFT. 
To this end, we have applied a small magnetic field of 
H = 0.005 eV, checked (for some temperatures) that 
this is still in the linear M vs. H regime (which fur¬ 
ther confirms the paramagnetic phase) and calculated 
x((f = 0) = M/H at this H. This way all vertex correc¬ 
tions are included; and this quantity allows to determine 
whether there is a second order phase transition towards 
a ferromagnetic phase or not. Prospectively competing 
phases with a different wave vector q are however not 
accessible this way. 


The full (f-dependent susceptibility x( 9 )*w) could in 
principle be obtained by solving the Bethe-Salpeter equa¬ 
tion. Unfortunately, this is computationally too de¬ 
manding for five orbitals at low temperatures. For 
the same reason the local susceptibility Xioc('r) could 
only be calculated reliably down to /3 = 30eV“^. 
But to gain at least some insight whether ferromag¬ 
netism or magnetic phases with other if-vectors pre¬ 
vail, we study the bare bubble susceptibility x'^ (?) = 

0) = ^)G'nmcr(*'4, fc -(- ()), 

which does not include vertex corrections. The result 
shown in Fig. 8 indicates that g = 0 is the leading 
instability.^^ Thus, in the following we will focus on 

x{q = 0). 
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FIG. 8: Susceptibility x°(' 7 ,= 0) vs. qx and qy (at qz = 
0) calculated from the DMFT d-electron Green functions at 
l3 = 100 eV~^. The maximum at q — 0 indicates that without 
vertex corrections ferromagnetism is the leading instability. 


FIG. 9: Ferromagnetic susceptibility x{<l = 0) as a function of 
temperature. The inset shows the temperature dependence of 
the corresponding local quantities: the local magnetic suscep¬ 
tibility (xioc(*cu = 0), diamonds) and the instantaneous one ( 
Xioc('r = 0), circles), from which the local magnetic moment 
can be estimated (oc ^xioc). 


The temperature dependence of the susceptibility 
xiq = 0) = M/H including vertex corrections is shown 
in Fig. 9. Upon decreasing temperature, we first notice 
an increase of the susceptibility. However, below 400K, 
the susceptibility decreases again. This clear trend of a 
reduction of the susceptibility by decreasing T makes the 
onset of a ferromagnetic order at lower temperatures ex¬ 
tremely unlikely. We note that a marked low-T reduction 
of x(q = 0) has been also reported experimentally^® and 
theoretically®® in the iron-pnictide compound LaFaAsO. 
There, this behavior of x(g = 0) coexists with an oppo¬ 
site (increasing) trend of the local magnetic susceptibility 
Xioc-^® Hence, the unusual low-T reduction of x(g = 0) 
has been attributed to specific features of the one-particle 
spectral function of LaFeAsO, displaying significant tem¬ 
perature variations near the Fermi level. By performing 
the same analysis for FeAl we find, however, that the low- 
T behavior of yioc (inset of Fig. 9) and xi? = 0) (main 
panel) is qualitatively very similar: both show a visible 
reduction for T < 400K. In the very same temperature 
interval, a slight reduction of the instantaneous local mo¬ 
ment (xioc(t = 0), inset) is also found, which is a typical 
behavior in the Fermi liquid regime, as described by the 
DMFT. 

Hence, in FeAl, the role played by emerging low-energy 
structures of the spectral function appears to be less im¬ 
portant than in LaFeAsO. Rather, the trend of x(g = 0) 
in FeAl may simply reflect the corresponding low-T re¬ 
duction of the local magnetic moment (oc y'xioc), hr par¬ 
ticular of the screened one. The latter can be ascribed to 
the enhanced metallic coherence of the low-temperature 
region, which is a general effect of local correlations in 
the Fermi-liquid regime. 


V. CONCLUSION 

In the present work, we studied FeAl by means of 
DFT-t-DMFT. In particular, we investigated the effects 
of electronic correlations on the magnetic properties since 
standard spin-polarized DFT calculations yield a ferro¬ 
magnetic ground state whereas in experiment it is a para- 
magnet. From DFT, we constructed a 9-band Wannier 
Hamiltonian with four A1 sp® orbitals and five Fe d or¬ 
bitals. For the latter we include a local 5'f/(2)-symmetric 
Kanamori interaction with U = 3.36 eV, U' = 2.36 eV 
and J = 0.71 eV as obtained from cRPA and solve the 
many-body problem by DMFT(CT-QMC). 

On the one particle level, the self energy and Green 
function suggest rather weak electronic correlations with 
a quasiparticle renormalization of only Z = 0.75 and no 
evidence of pronounced upper and lower Hubbard bands. 
In DFT-I-DMFT we calculate from the spin-spin corre¬ 
lation function an equal-time local magnetic moment of 
1.6 pb which is twice as large as the magnetic moment 
in spin-polarized DFT. It is also much larger than the 
bubble contribution which demonstrates that electronic 
correlations are pivotal for the two-particle quantities in 
FeAl. 

Even more importantly, our results show that the mo¬ 
ment is fluctuating in time and screened on the fs time 
scale. This explains why we also do not find long-range 
ferromagnetic order. According to our DFT-I-DMFT 
study, FeAl is paramagnetic with a maximum in the fer¬ 
romagnetic susceptibility around room temperature and 
no tendency towards long range magnetic order in the 
temperature range studied. 

Previously, it has been proposed that disorder and a 
spin-glass behavior might explain the missing ferromag- 
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netic moment in experiment. Our results show that if 
temporal fluctuations are taken into account the moment 
is actually screened on short time scales. There is hence, 
even for a perfect lattice, neither ferromagnetism nor a 
local moment on longer time scales. If the magnetic mo¬ 
ments were constant in time and spatially disordered, 
Mossbauer experiments, which probe the local magnetic 
moment, are in principle able to reveal it. The fluctu¬ 
ating local moment on the fs time scales, can however 
not be observed in Mossbauer spectroscopy which can¬ 
not resolve such short time scales. Hence, the Mossbauer 
experiments,^ which show no local magnetic moment for 
stoichiometric FeAl, seem to better agree with a magnetic 
moment fluctuating in time as we find in DFT-I-DMFT 
than with a moment fluctuating in space. 
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